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Abstract 

Integral membrane proteins deform the surrounding bilayer creating long-ranged forces that 
influence distant proteins. These forces can be attractive or repulsive, depending on the proteins' 
shape, height, contact angle with the bilayer, as well as the local membrane curvature. Although 
interaction energies are not pairwise additive, for sufficiently low protein density, thermodynamic 
properties depend only upon pair interactions. Here, we compute pair interaction potentials and en- 
tropic contributions to the two-dimensional osmotic pressure of a collection of noncircular proteins. 
In contrast to direct short-ranged interactions such as van der Waal's, hydrophobic, or electrostatic 
interactions, both local membrane Gaussian curvature and protein ellipticity can induce attrac- 
tions between two proteins at distances of up to ten times their typical radii. For flat membranes, 
bending rigidities of ~ SOksT, and moderate ellipticities, we find thermally averaged attractive 
interactions of order ~ IUbT. These interactions may play an important role in the intermediate 
stages of protein aggregation. 
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>1 Introduction 

^ |Membrane proteins interact directly via screened electrostatic, van der Waal's, and hydrophobic forces. 
■ - - 'These are short ranged, operating typically over distances of less than a nanometer. Proteins can also 
interact indirectly via the bilayer in which they are dissolved. In particular, a protein that is "geo- 
metrically mismatched" to the bilayer will induce deformations that affect neighboring proteins. These 
"solvent induced forces" (the membrane lipids being the solvent) are generated by bending deformations 
of the bilayer and typically act over many protein diameters. 

Membrane associated proteins can aggregate due to bilayer bending mediated interactions. For 
example, aquaporin AQPl and CD59 aggregate to tips of pipette-drawn tubules (|Cho et al., 1999 



pischer fc Mohandas, 1996| ). Previous studies using a continuum approximation for the intervening 
bilayer membrane, have treated protein-protein interactions and found an repulsion between two 
identical inclusions ( [Goulian, Bruinsma fc Pincus, 199^ ; |Kim, et al., 199B| |Park fc Lubensky, 1996 
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Dommersnes, Fournier fc Galatola, 1998|) . Goulian et a/.( |Goulian, Bruinsma fc Pincus, 1993|) also find 



a weak attractive {—ksT/r^) interaction arising from Casimir forces resulting from suppressed thermo- 
dynamic fluctuations of the intervening membrane. Here, we study in detail a direct mechanical origin 
for protein-protein attractive interactions. Although bending induced forces between multiple inclusions 
are not pairwise additive, ( [Kim, et al., 199^ ; [Kim, et al., 1999| ), in this paper, we restrict ourselves to 
low protein densities where only pairwise interactions are relevant. We find that the interplay between 
protein noncircularity ( [Kim, et al., 199"9| ) and background Gaussian curvature curvature dramatically 
affect protein-protein attractions and thermodynamics. 

Many membrane proteins are noncircular in the plane of the membrane, including adsorbed polypep- 
tides such as MARCKS (|Myat et al. 19971) , and bacteriorhodopsin ( [Luecke et al. 1999[) , which consists of 
seven transmembrane helices arranged in an elliptical configuration. Small domains, dimers, or droplets 
of molecules such as cholesterol or specific lipids can themselves behave effectively as membrane in- 
clusions. Droplets need not be rigid to induce attractive forces among themselves; fluctuations in the 
droplet shape itself may lead to an effective attraction. 

In Section 2 we briefly review the mechanical theory of inclusion-induced bilayer bending ( [Helfrich^ 
1973 |; |Kim, et al., 1998[ ; [Netz fc Pincus, 1995[ ). The lipid membrane is approximated by a thin plate 
that resists out-of-plane bending. Inclusions such as integral membrane proteins, or surface adsorbed 
molecules, impose boundary conditions along the contact line between the membrane and the protein. 
Using elastic plate theory to describe the membrane deformations, we derive the energy for two identical 
inclusions as a function of their relative position within the membrane surface. 

In Section 3, we show that the rotational and translational time scales can be separated, so that 
we can thermally average out the fast rotational degrees of freedom. The resulting effective potential 
between two proteins is attractive provided the inclusions are sufficiently non-circular. We use the 
effective potential to compute the second virial coefficient and show how the attractive interactions 
affect the two-dimensional protein osmotic pressure. Finally, we discuss biological processes where 
membrane induced long-ranged protein-protein attractions may play an intermediate role, and propose 
possible measurements. 



2 Membrane inclusions and height deformation 

Small membrane deformations (on the scale of the lipid or protein molecules) can be accurately modeled 
using standard plate theory ( [Landau fc Lifshitz, 1985[ ; Pelf rich, 1973[) 



E[H{S), K{S)] =2b^ dSH\S) + hg j dSK{S), 



where H{S) and K{S) are the local mean and Gaussian curvatures, and b and bg are their associated 
elastic moduli. We have assumed a symmetric bilayer and a vanishing spontaneous mean curvature in the 
absence of the membrane-deforming proteins of interest. For uniform bg, the Gaussian contribution (the 
second integral in Eq. |l|), when integrated over the entire surface yields a constant that is independent 
of the relative configurations of the embedded proteins ( [Kim, et al., 1995| [Struik 1994[ ). Thus, the 



Gaussian energy term can be ignored when considering protein-protein interaction energies. 

In an expansion of the free energy about that of a fiat interface, H{S) ~ (l/2)V^/i(x, y), where is 
the two-dimensional, in-plane Laplacian, and h{x,y) is a small, slowly varying height deformation from 
the fiat state (cf. Fig 1). Minimizing E[h{S)] with respect to h{x, y € S) yields the biharmonic equation 



V^h{x,y) = 2V^H{r) = 0. 



(2) 



First consider membrane deformations about an isolated, circularly symmetric inclusion of radius a. 
If the bilayer midplane contact^ the protein perimeter C (see Fig. |l|) at an angle 7, the appropriate 
solution to Eq. ^ is h{r) = — 7ln(r/a) for r > a. We have excluded terms in h{r) of the form 
In r, r^, const, because they are unbounded in energy (Eq. |l]), or do not satisfy the contact angle 
boundary condition at r = a. Since V^ln(r/a) = 2H{r) = for r > a, there is no mean curvature 
bending energy (proportional to b) residing in the bilayer. 

Now consider cases where more than one inclusion are present, or where the contact angles, heights of 
contact, or the shapes of the membrane associated proteins are noncircular. Three types of noncircularity 
can arise. The inclusion itself may be noncircular {e.g. elliptical), the height of contact of the bilayer 
midplane to the inclusion may vary along the perimeter C of the protein, and the contact angle itself 
may vary along C. These noncircular boundary effects arise from the detailed microscopic nature of 
the protein and its interaction with the lipid molecules. When more than one protein is present, the 
deformations surrounding each protein are not circularly symmetric. A nonvanishing mean curvature, 
H{r), that gives bounded bending energies can be represented by a multipole expansion, 

00 

if (r, e) = ^ r~" (a„ cos nO + sin n9) , (3) 

n=2 

where (r, 6) is the radial and angular coordinate about an arbitrary origin. Upon substitution of Eq. ^ 
into Eq. 0, we find the bending energy E ^ b Yl^=2^'^'n + ^n)- To determine a„, 6„, we solve V^/i(r, 6) = 
H{r, 9) and impose boundary conditions (see Appendix A) on /i(r, 9) at C. In the limit of small 
noncircularity or low protein concentrations, the largest nondivergent terms are associated with n = 2. 
Wiggly inclusion cross-sections or highly oscillating boundary conditions only weakly affect membrane 
bending-mediated protein-protein interactions via n > 2 terms. We derive the full multibody, interaction 
energy in Appendix A. The two-body interaction energy measured in units of ksT is 
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E{R,9u92;A,K,,n)= ^ + K,--e~'''^ + ^ + K, - -e-'^'^ . (4) 
The dimensionless separation distance R, protein ellipticity A, and background curvature Kb are given 

by 



where B = nb/kBT is the dimensionless bending stiffness, and e ~ 0{e) quantifies the noncircular 
nature of each inclusion (see Appendix A). The angle Q is measured between the line joining the protein 
centers and the principle axis of curvature defined by the background Gaussian curvature (see Fig. |^). 
The angles 9i , ^2 are measured between the principle axes of proteins 1 and 2 and the same principle 
axis. The quantity Kb measures the local, externally induced (via other distant proteins or external 
bending forces) background curvature in this principle axis direction. We show in Appendix A that the 
dominant effect of distant proteins is to induce mean curvature deformations that decay as but 

^The contact angle 7 incorporates the details of the molecular interactions between the included/ adsorbed protein with 
the lipid molecules. A molecular dynamics simulation would in principle provide the basis for a quantitative estimate of 
7, but is beyond the scope of this paper. We will simply assume that 7 is a phenomenological parameter determined by 
the local chemistry, in complete analogy with the standard liquid-gas-solid contact angle. 



constant negative Gaussian curvatures. The local curvature Kb arises only from deformations that are 
of zero mean curvature. In what follows, our statistical thermodynamic analyses will be applied to the 
pair interaction energy given by Eq. ^ with the convention A,Kh > 0. 



3 Rotationally averaged interactions 

Proteins that are not attached to the cytoskeleton are free to rotate and diffuse within in the membrane. 
The interaction potential between two membrane-deforming inclusions is a complicated, nonseparable 
function of their relative major axis angles and separation distance (cf. Eq. Although the energy is 
a function of the specific separations and angles between two membrane associated proteins, rotational 
degrees of freedom are sampled faster than the translational of freedom. This can be shown by the 
following argument. 

A small solvent molecule in solution has a rotational correlation time of the order Trot S 1 ^s, while 
its translational diffusion constant is Dt^uns ~ 10~^cm^/s. Therefore, in the time r^ot it takes for a small 
solvent molecule to lose rotational correlation, it would have translated 

5r ~ \/TrotDtrans ~ O.lum. (6) 

For membrane constituents, such as bilayer lipid molecules. Trot ~ 1 — 5ns, and Dtrans ~ lO^^cm^/s, 
where Trot corresponds to rotation about the molecular axis parallel to the normal vector of the membrane 
([Marsh, 1990|) . As with small molecules in bulk solution, membrane-bound lipid molecules also move 



6R ~ O.lnm during a rotational correlation time. For larger membrane inclusions such as proteins, both 
rotation and translational diffusion are slower. If A is the relative size of the protein radius with respect 
to the effective lipid radius, protein rotational correlation times increase by ~ A'^ while Dtrans decreases 
with a. Membrane proteins that are not too large, A < 10 say, diffuse 6r ~ Inm during the time 
over which it has lost rotational correlation. Therefore, in the time it takes for a typical inclusion to 
rotate about its axis, it has diffused less than its own size, typically >lnm. This estimate is consistent 
with fluorescence measurements that find Trot ~ 0.1 — 1ms ( [Yamada et al. 1999] ). Rotational time 



scales for larger proteins may not be much faster than translational motions, therefore, our subsequent 
model is most appropriate for small, unhindered membrane proteins. We must eventually verify that 
the protein-protein separation r of interest is greater than the typical diffusion distance 6r. 

The time scale separation can be implemented by statistically averaging over the principle axis angles 
of the two inclusions while keeping the distance R and angle Q between them fixed. Weighting the exact 
two particle energy over its own Boltzmann weight, 

E,ff{R; A, Kb, n) = / E{R; A, K^O,, ^2)e-^(^''^''^^^'^-^)(i^i(i^2, (7) 
where the rotational partition function 

Z= (8) 
Jo 

Upon substitution of Eq. | into Eqs. |^ and H, and performing the integration (see Appendix B), 

S,^(i?;ir„n) = ^ + ^ -2^111 (9) 

where 



The effective interaction potential between two inclusions is defined by tlie difference between tlie mem- 
brane bending energies of two inclusions separated at distance R and at infinite separation, 



U,ff{R- A, Kh, n) = E,ff{R- A, K,, Q) - E,ff{oo) 



2e 2e/i(0 



(11) 
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For fixed ellipticity A, the set of parameters Kh, and R that gives rise to a minimum at R* < oo, if it 
exists, is implicitly determined by 



dR 



0. (12) 



For sufficiently small R, Ueg — 2/-R^, as in the circular protein case. 
3.1 Zero background curvature 

First consider the case of two isolated proteins embedded in a fiat membrane. In the absence of external 
mechanical forces that impose background membrane deformations, and with other inclusions sufficiently 
far away, Hb = Kb = 0, and ^ = |A|/i?^. The effective potential (Eq. |lT]) becomes 



t/.,(«;A,.-.^0).i,-(^)^. (13) 

Without background curvature {Kf, = 0), there are no defining principle axes, and Ueff is independent 
of how the angle of the segment joining the inclusion centers is aligned. Clearly, an effective attractive 
interaction can arise for A/R^ > 1, when h{A/R^)/Io{A/R^) ~ 1, and Ueff{R]A,Kb = 0) ~ l/R^ - 
A/R^. Although the interaction (Eq. |) yields both repulsive as well as attractive forces, the Boltzmann 
thermal average in Eq. |^ favors the lower energy configurations of 61,62- Hence the pair of inclusions 
spends more time in attractive configurations, resulting in a residual attraction in Ueff{R)- In the Kb = Q 
limit, the large R behavior of Eq. |TB| is 

U.ff{R) = '^-^ + 0{R-'). (14) 

Since the potential becomes repulsive at short distances, an effective ellipticity A > A* = \/2 is necessary 
for the existence of a minimum in Uejj{R). 

Figure ^ shows the 6'— independent effective interaction potential as a function of R for various A. 
As A is increased from A* = v^, the minimum radius R* determined by Eq. 0, decreases rapidly from 
R* ~ 00. The A > A* dependence of R* is plotted in Figure ^d. Also shown are the corresponding 
magnitudes of the global minima of Ueff{R; A, Kb = 0) as a function of A. 

Figures 0a,b show that appreciable attractive wells can persist at distances i? ~ 1. For example, 
the minimum determined by the set of parameters A ?s 2, i?* ?a 0.9 and \Ueff{R*)\ ~ 2(A;bT) can 
arise for 7 ~ 1, e/a ~ 0.25 and b ~ SOksT. The separation corresponding to the minimum energy in 
this case is r*/a ~ 9, or nine times the inclusion radius. Our initial assumptions using measured and 
estimated rotational/translational diffusion constants for typical membrane proteins are validated since 
5r/a~l — 2<C r*/a. We conclude that thermally averaged noncircular membrane deformations can 
induce long-ranged attractive interactions of at least ~ 2kBT up to distances ~ 10a. 



3.2 Effect of local Gaussian curvature, Hb = 0, 



A local background curvature may arise due to a nonuniform distribution of distant membrane proteins 
or an externally imposed deformation. For example, in the experiments of ( piio et al., 1999| ; pischer S 



[Mohandas, 1996|) , a cell is manipulated by a micropipette. A lipid neck is drawn into the pipette and 
curvature is externally imposed. Regions near the base of the neck will have a large negative Gaussian 
curvature. Similarly, membrane fusion and fission processes in endo/exocytosis involves intermediate 
shapes with constricted necks containing Gaussian curvature. These regions may be "externally" im- 
posed by proteins involved in vesiculation {e.g. dynamin or motor proteins). The Gaussian curvature in 



this case may also result from lipid structural or composition changes ([Schmidt et al.., 1999|) . Therefore, 
curvature can couple to membrane protein or lipid shapes. 

The Gaussian curvature of the membrane between the two proteins establishes local axes of principle 
curvature such that ad'^_^h{xi,X2) = —ad'^^h{xi,X2) = rjf, oc Kb > 0. Since we assume Hf, = 0, the back- 
ground deformation between the two proteins will resemble a saddle with principle curvatures of equal 
magnitudes (cf. Fig. The rotationally averaged (over 6*1, 6*2) effective interaction Ueff{R;A,Kb,Q) 
will generate attractions at specific orientation angles Q even if A < A*. This can be most easily seen 
by expanding equation |ll| (the rotationally averaged interaction Ueff{R; A, Kj,, fl)) in powers of for 
large R: 

U,ff{R ^ 00; A, Kb, fi) ~ ^ cos 2fi + ^ + O , (15) 

where explicit forms for A2, A4 are given in Appendix A. The appearance of 7^ when Kb > 
immediately generates a minimum, however small. Even when ellipticity vanishes (A = 0), oc 
Kb cos2Q < for appropriate Q. 

The physical origin of attractions in the presence of background curvature can be readily seen by 
considering Figure |^. With our convention 7 > 0, circular proteins situated at low regions of the saddle 
{Q ~ 7r/2) develop attractive interactions, while those with i7 ~ always repel. Recall from previous 
studies that two circular protein repel with a i?"^ potential ( Poulian, Bruinsma fc Pincus, 1993[ ; [Kim, et 



M-1 1998 ; Park fc Lubensky, 1996 ; Dommersnes, Fournier fc Galatola, 1998 ). This is a direct consequence 



of placing a second protein in the Gaussian curvature of the first one. When the background curvature of 
the membrane in the region between two proteins augments the individual Gaussian curvatures around 
the first protein (near = 0), the R~^ repulsion is also enhanced. Conversely, if the background 
curvature mitigates the saddle induced by an individual inclusion (near Q = vr/2), the other inclusion 
sees not only a diminished repulsion, but a mutual attraction at large enough distances. This is because 
the individual Gaussian curvature around a protein (arising from h{r) ^ — 7ln(r/a)) decays as 1/r^ and 
eventually becomes smaller than the imposed constant background Gaussian curvature associated with 
Kb. Attractive effects of the background curvature eventually manifest themselves when Q ~ 7r/2. 

Figure ^ shows the effects of a small amount of local background curvature on the effective interac- 
tion potential. For small ellipticity A -C A*, minima appear for large enough angles Q (approximately 
for Q > it/ 4:). For similar background curvatures, but much larger ellipticities, the potential develops a 
repulsive barrier before becoming attractive for certain Q. This signals that ^4 < for large enough A 
and is depicted in Fig. for A = 2.5. In the limit of small Kb, < when 

A > A* + ^ + ^(3 + sin^ 2Q)j + O (K^) . (16) 

There is yet an additional, qualitatively different feature of Uejf{R', A, Kb, f2) when both A and Kb are 
large. Although typical values of Kb (see Eq. H) in biological settings is Kb -C 1, we find that large 
values of Kb and A give rise to double minima in the interaction potential, especially near Q ~ 7r/2. 



Figure ^ shows double minima for Q = Ttt/IG, 7r/2. Additional higher order coefficients such as Aq/R^, 
etc. are required to quantitatively describe multiple minima. The two minima are a consequence of 
the two independent physical effects that prefer energy minima; local Gaussian curvature associated 
with Kh and effective ellipticity A. Typically, the weaker, longer-ranged minimum is predominantly the 
signature of a large Kh, while the deeper, shorter-ranged minimum (such as that shown in Fig. |^ and 
|b) is a feature of ellipticity A > A*. Saddles of order Kb > 1 correspond to principle radii of curvature 
on the order of ~ 10 times the protein size a, and are thus regions of extreme Gaussian deformations. 
Regions of such warp may would only exist as transient, small systems such as fusion necks. Henceforth, 
we will restrict ourselves to Kf, small enough to only induce one minimum. 

Angles Q which yield attractive interactions can be estimated by considering ^2,^4. Assuming 
A4 > 0, values of < give attractive interactions when — 7r/4 < ^7 < 7r/4. When A2 > 0, proteins 
with orientation 7r/4 < \Q\ < Sir /A will experience attractive forces. However, these conditions are 
modified if A4 < 0, when some angles within —ir/A < Q < n/A can yield attraction even if A2 > 0. This 
case corresponds to Fig. ^ where a repulsive barrier ai R> R* arises. A minimum can still arise even 
at angles where A2 cos 2^2 > due to the —R~^ behavior. The matching to repulsive behavior at smaller 
R requires consideration of +R~^ terms. 

The top panel of Figure |^ shows the radius corresponding to the only minimum of the effective 
potential Uejj as a function of Kb, for A = 0.5, 2 and 4. Both east-west and north-south configurations 
are shown, with intermediate angles Q interpolating between the curves. For small ellipticity, the local 
principle curvature Kb is the predominant source of attraction at larger distances, shown by the thick 
dashed curve. Increasing Kf, destabilizes the effective energy minima near i7 = 0. Above a certain 
background Gaussian curvature intensity, the effective potential minimum evaporates to R* —>■ 00 for 
proteins situated at = (solid curves), and the attraction is washed out. For small Kf,, the two effects, 
ellipticity and background Gaussian curvature, complement each other near Q = 7i/2 in reinforcing an 
energy minimum. Consistent with Fig. |a for A > V2, R* in Fig. | (thick curves) is smaller for larger 
A. The bottom panel plots the corresponding minimum energies. 

The fi-dependence of R* and the minimum energy is shown in Figure ^. As expected, or large 
A ^ a/2, both R* and Ueff{R*,^) are fairly insensitive to il. When A is small, the energy minima 
and their associated radii R*, caused predominantly by Kf,, are very sensitive to orientation Q. These 
behaviors are consistent with the energy profiles shown in Fig. ^3. In fact, for small enough A, the 
minima near ^ are annihilated, independent of Kb. Thus, we see a qualitative difference between 
attractive potentials generated by intrinsic ellipticity and background Gaussian curvature. 



4 The second virial coefficient 

We now consider the influence of the effective protein-protein attractions on a low density ensemble of 
inclusions. By analogy with the molecular origins of the osmotic second virial coefficients of proteins 
in solution ([Neal, Asthagiri fc Lenhoff, 19981) , we will consider the bending energy contributions to 



the second virial coefficient for a two-dimensional protein equation of state. The membrane mediated 
interactions however, are much longer-ranged than those in solution ( [Neal, Asthagiri fc Lenhoff, 1998| ). 
Consider the thermodynamic limit and times long enough such that 

r > > Trot, (17) 

''-^ trans 

where Dtrans is the protein translational diffusion constant. On the time scale r, the inclusions are 
relatively free to diffuse about the bilayer. They interact among themselves via the rotationally averaged 
potential Uejj that manifests itself on time scales > Trot- For very low protein densities (large protein 



separation £), the two dimensional protein osmotic pressure will be nearly that of an ideal gas, analogous 
to a low density gaseous phase surfactant monolayer at the air water interface. Finite protein size a, 
and longer-ranged elastically-coupled interactions will give nonideal gas properties. The first correction 
to ideahty in the equation of state is given by the second virial coefficient (|McQuarrie, 19761) : 



^ r + + O (BsT^) 



(19) 



where T is the surface concentration of protein and B2 is computed using the formula 

-| /*27r POO 

2^ Jo Jo 

1 p27r POO 

= -- / (e-^<=i^(^^^'^'"^)+^-(^'^'')-l)i?di?cifi, 

2 Jo Jo 

The second virial B2 represents the small fraction of pairwise interacting proteins. Equations and |TP| 
are nondimensionalized such that the surface density F -C 1 is measured by the number of proteins in 
area (see Eq. ^ and the protein osmotic pressure 11 is measured in units of ksT/Rl. Equation |1^ is 
exact and does not require the separation of rotational and translational diffusion times needed for the 
derivation of Uejf{R] A,Kb,Q). Here, we do not consider how integrating out the rotational degrees of 
freedom affect the fixed translational degree of freedom. Instead, we are considering times long enough 
for equilibration of both degrees of freedom, and their combined contribution to the equation of state 
via B2. 

The physical origin and value of Kh used in Eq. |19| is as follows. The local curvature felt by the 
interacting pair represents an interaction between this pair and some other distant proteins. However, 
the virial equation of state (Eq. p^) is a systematic expansion in surface density expanded about an 
ideal, noninteracting ensemble. Since membrane bending mediated interactions are not pairwise additive 
(|Kim, et al., 1998]) , one might be tempted to assume that the presence of other proteins would modify 



the interaction energy E used in the expression for i?2. However, these more complicated interactions 
would depend upon the concentration of the other background proteins, and would generate terms 
higher order in F. In other words, we start at densities so low that the protein ensemble is completely 
noninteracting. As the density is slightly increased, a pair of protein molecules occasionally interact 
and perhaps form dimers, with each pair ignorant of any other protein. At this still rather low density, 
the probability three or more proteins approach each other is negligible. When the density is further 
increased, one needs to consider the higher order virial terms. Therefore, to second order in F, the 
deviation of the equation of state from ideality is completely determined by the two-body interaction 
E{R,9i,92;A,Kb,Q) and is independent of nonpairwise effects ( pVIcQuarrie, 1976| ). Note however, that 
the two-body interaction will depend on the Kh associated with externally forced, zero mean curvature 
membrane deformations. Therefore, for the expansion Eq. |T8|to be consistent, the value of Kf, = Kh to 
be used in Eq. |l^ is that owing solely to external force-generated Gaussian curvatures, independent of 
the protein density. 

Figure ^ shows the numerically computed second virial coefficient as a function of inclusion ellipticity 
for various Kh. As expected for small Kh, the virial coefficient becomes increasingly negative as the 
ellipticity increases. The value for circular inclusions B2{A = 0, Kh = 0) = vr'^/^/y^ corresponds to 
purely repulsive disks with mutual interaction U{R) = 2/R'^. At ellipticity A 1.69, -82(1-69, A'b = 
0) ~ corresponding to a protein solution that is ideal to second order in surface density. Although 
when A ~ 1.69 > A* = ^2, Ueg has an attractive minimum, its effects are negated by the repulsive R 
part of the interaction such that the overall, effective contribution to B2 vanishes. For A > 1.69, the 
effective attraction between membrane proteins begins to manifest itself. The second virial is modified 



by externally imposed Gaussian curvature. Recall that when Kb ^ 0, certain angles lead to attractive 
interactions, even for small A < A*. Since we are now thermodynamically averaging over protein 
positions and f2 in addition to ^1,6*2, the inclusions will spend more time at attractive, lower energy 
angles f2, hence lowering i?2. Consistent with Fig. ^, larger values of for A > A* lead to stronger 
repulsions at small fi, which average into i?2, making it less negative. 

The dependence of Bi on is indicated in Figure l^^- In the absence of ellipticity, B2 is given by 
the integral 

52(A = 0;i?b) = -TT / e~^'^'h{AKi,/R^) - 1 RdR. (20) 

JQ L -I 

For A > 0, i?2, found from numerical integration of the full expression Eq. are also shown in Fig. 

For Kb = 0, increasing ellipticity decreases inclusion repulsions and i?2- As in Fig. 0a, large Kb and 
A tend to increase -82- 

Equation | was used in Eq. to compute the curves shown in Figures thus, the protein- 
protein interaction was assumed to consist of contributions only from membrane bending. The hard 
core, excluded area of each protein, ~ vra^, can be included by modifying Ueff{R) by setting Ueff{R < 
a/Ro) = 00. Although we expect this additional repulsive term to further reduce the effective sampling 
area of the inclusions, and increase the second virial coefficient, we find that for all reasonable values of 
Rq, B2 does not change noticeably from those shown in Figs. |^. The hard core part of the potential, 
due to e.g. close-ranged van der Waals repulsion, is not statistically sampled by the inclusions since the 
membrane bending induced interactions (~ 1/-R^) already keeps them far apart. 

Since nonpairwise interactions manifest themselves only at third and higher order in F, we can 
estimate their importance by comparing B2T^ with -BsF^. For nonpairwise interactions to be thermo- 
dynamically relevant it is necessary but not sufficient that the surface density 



r > 



Bo 



(21) 



i?3 

Although multibody interactions may be important microscopically, their effects on the low density 
equation of state, cannot be resolved. Even if the density is high enough for B^V^ to be measurable, the 
value of i?3 is found via a four-dimensional integral over configurations of three membrane proteins. All 
orientations and distances will be averaged and all components of their interactions, repulsive, attractive, 
pairwise, and nonpairwise will be included. In other words, one cannot uniquely determine the potential 
U from a measurement of B^. 



5 Discussion and Conclusions 

Proteins beyond the range of screened electrostatic, or van der Waals molecular forces can exert forces 
on one another by virtue of the deformation they impose on the lipid bilayer. These interactions can be 
attractive if the proteins have a noncircular cross sectional shape or if the local membrane deformation 
is saddle shaped (negative Gaussian curvature). For bending rigidities h ~ SOksT, and protein shape 
ellipticities e/a ~ 0.25, we find attractive interactions of a few fc^T acting at a range of ~ 5 protein 
diameters. Thus proteins of radii ~ 5nm can interact at distances of ~ 50 nm, much further than any 
interaction between similar molecules in solution. On a fiat membrane {Hb = Kb = 0), an effective 
ellipticity e > {2kBT /nhY^'^ is necessary for a potential minimum to emerge between a pair of proteins. 

Although we have presented the model in terms of integral membrane proteins, noncircular pe- 
ripheral proteins, or lipid/cholesterol/proteins aggregates can also induce local membrane bending. 
Elastically-coupled interactions between peripheral membrane proteins can mediate dissociation and 



binding. Dimerization of noncircular peripheral proteins lowers their absolute energy below that of sep- 
arated ones, and so they are less likely to dissociate from the membrane. Similarly, in one dimension, 
adsorbed proteins can bend DNA and affect the binding of a second nearby protein. The effective in- 
teractions in that case also depends upon the protein orientation ([Diamant fc Andelman 1999| ; [Rudnick 
& Bru jinsma, 1999| ). Moreover, proteins and protein aggregates need not be rigid, as we have assumed 
here. Noncircular distortions of a lipid or protein domains can fluctuate in such a way as to yield 
domain-domain interactions of exactly the same form considered here. 

Elastically mediated attractions can also manifest themselves in the aggregation of circular proteins. 
Once circular proteins overcome short-ranged repulsions and dimerize, barriers to further aggregation of 



these elliptical dimers are reduced by dimer-dimer attractions described by Eq. ([TT]) . If the inclusions 
are themselves dimers or higher aggregates that persist on the time scale of rotation, bending mediated 
attraction would enhance further aggregation. 

We have only considered the mechanical energies of the intervening lipid bilayer. However, an 
ensemble of membrane-bound proteins or a mixture of lipids can manifest itself through other forces. 



For example, the presence of charged membrane components can induce bending (|Chou, Jaric fc Siggia. 



1997|2and initiate endo/exocytosis and organelle trafficking. Protein-protein interactions arising from 
screened electrostatic forces operate at much shorter distances than those of the bending induced forces. 
Therefore, by spatially organizing charged membrane components, elastic interactions may also play an 
indirect role in large scale electrostatically induced membrane deformation. 

We also considered an ensemble of surface proteins elastically coupled by membrane deformation 
and computed the deviation of its equation of state from that of an ideal solute. Although membrane- 
mediated protein-protein interactions are nonpairwise additive (|Kim, et al., 1998|) , only the two-particle 
interaction is relevant for sparsely distributed proteins. On a flat membrane the second virial coefficient 
B2 < when e > Q.9b\JkBT /b (A 1.69). At this ellipticity, the elastically induced repulsive 
interactions just compensate for the rotationally averaged attractions. This dependence on h/ksT 
suggests that the cell can regulate protein-protein interactions by varying the lipid composition and 
hence the bending rigidity of the bilayer. Thus, the formation of cholesterol rafts and lipid domains may 
have an indirect role in mediating long-ranged surface protein aggregation and activity. 

Finally, we propose possible experiments in artificial membrane systems where the surface density 
can be made small enough for a virial expansion to be valid. Although the two-dimensional osmotic pres- 
sure would be difficult to measure accurately, measurements of the association time between dimerized 
proteins are feasible. Measurements have been made of the lifetimes of gramicidin A channels composed 
of dimers of barrels in opposite bilayer leaflets as a function of bilayer thickness ([Kolb fc Bamberg 1977| ; 
[Elliot et al. 1983]) . Measurements of dimer lifetimes as a function of lipid tail length may reveal the de- 
pendence of the attractive interactions outlined in this paper. In fact, since the second virial coefficient 
measures the time-averaged fraction of proteins in dimers at low density, their lifetimes are proportional 
to B2 for attracting proteins. Moreover, these association lifetimes can be measured in the presence of 
externally imposed Gaussian deformations. Even though an imposed Gaussian curvature increases the 
interaction well depth at Vt ^ 7r/2, and destroys the attractions for proteins near f2 0, the overall 
statistical effect, is to enhance binding, as is evident from Figures |^. Therefore, we expect that dimer 
lifetimes can be enhanced for proteins residing in regions of large magnitudes of Gaussian curvature 
such as the base of extruded tubules. This may be instrumental in recruiting fusagens to the correct 
location for membrane budding. Dimer lifetimes potentially can be measured by fluorescence transfer 
of specifically designed hydrophilic moieties attached to membrane proteins. Nonpairwise interactions 
can only be probed directly by measuring lifetimes and aggregation rates of trimers. This would require 
statistical analyses of chemical or fluorescence activity among two differently tagged membrane proteins, 
or single molecule diffusion studies. 
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A Interaction energy among noncircular inclusions 

We consider the boundary conditions that the height, h{r, 6), must satisfy and the effects of noncircular 
proteins on the interaction energies ( [Kim, et al., 19991 ). Consider proteins with chemistry that changes 
the cross-sectional protein shape from circularity by an amount e. The concomitant changes in lipid 
contact height and angle are also assumed to be modified by 0{e). As shown in Fig. |l|, the protein 
perimeter, measured from the protein center is, to order 0{e), 

C-{a + ecos2{e-9i))n, (22) 

where n is the unit normal vector to the curve C projected onto the bilayer midplane, and £:cos2(6' — ^j) 
is a small, angle-dependent perturbation measuring the deviation from circularity of protein i. Upon 
expanding the general boundary conditions h{C) = 6h{6) and n ■ Vh{C) = —7 — S'y{6) to lowest order 
in e, we arrive at effective boundary conditions: 

h{a) ^ 6h{e - 6,) + 75 cos 2{e - Oi) + 0{e^) 

/ e ^ ^23) 

drh{a) ~ -7 (^1 + - cos 2(0 - 9i)j - 5-fi9 - 9i) + 0(8^) 

where we have for simplicity also assumed the variations 6h{C) and 57(C) to be also of order e. 

In the limit of small noncircularity or low protein concentrations, the dominant nondivergent con- 
tribution of H{r) to the energy is + The deformation h{r, 9) that satisfies V^/i(r, 9) = 2H{r, 9) 
and Eqs. ^ can be written in the form 

00 

h{r, 9) ~ —7 In ^-j + ^ {fn{r) cosn9 + gn{r) sin n6') , (24) 

^ n=2 

and determine 02, &2- When the proteins have intrinsic noncircularity (e 7^ 0), 02 + ^2 turns out to be the 
magnitude of the local Gaussian curvature (since Hb = 0), modified by additional 6'— dependent terms 
(|Kim, et al., 19991) . The local Gaussian curvature due to the other j far field proteins, in either case. 



is calculated using the leading order term h(r) ^ —7 In \r — fj\, which is simply a superposition of the 
longest-ranged Inr terms about each inclusion. The total bending energy E[H{r, 9)] for an ensemble of 
inclusions can be written in the complex form ( [Kim, et al., 1999 ), 



a2 



^ = -^^^E Ezr^T^-^e--^- • (25) 



where Zi = Xi + ii/i is the position of the i*^ protein in the complex plane, and 



i)(j + 2^--6j] (26) 



a/ \ a 



measures the effective ellipticity of the identical proteins. Now consider two relatively isolated, identical 
proteins i,j = 1,2. The effects of proteins far away are felt via a local Gaussian curvature emanating 
from these background proteins. Upon explicitly separating these contributions, the pair interaction 
energy becomes 



+ 



+ r/, - -e 



where 



Vb = a- 



dxj 



(27) 



(2^ 



is the curvature in the xi principle direction due to far-field background inclusions or externally induced 
deformations hf, ~ — 7ln|z — Zj\, j > 3. The mean curvature expanded about a noncircular protein 
(Eq. ^) results in a deformation h{r, 6) with terms proportional to cos 26, sin 26 ( |Kim, et al., 1998| ). 
These term carry zero mean curvature, but constant negative Gaussian curvature. From the expansion 
Eq. ^ the only mean curvature contribution can be seen to decay as r~^, which we neglect. A further 
contribution to the local saddle curvature, 7]^, felt by the two proteins can arise from externally applied 
mechanical forces that deform the bilayer in an appropriate way. The angles 6i , 62 are the angles of the 
principle axes of the inclusion shape (or the height or contact angle functions Sh, ^7) measured from 
the principle background curvature axis xi. The angle Q measures the angle between the principle 
background curvature axis and the segment joining the centers of the two inclusions. Upon rescaling 
according to Eq. ^, we arrive at the energy given in Equation ^. 



B Rotational averaging 

The integrals 

/ E{R,6i,62;Kb,n)e-'^d6id62 and Z^/^ = / e-^d6id62 (29) 
Jo Jo 

used to compute the rotationally averaged, effective protein-protein interaction involve integration of 



y gacos20+/3sin2e^^ and ^(a COS 2^ + /? siu 2^)e"™^'^+^""2^(i^. 



(30) 



The first integral in Eq. [3^ can be computed in closed form by substituting the exponents with their 
Bessel function expansions 



„a cos 29 



/o(a) + 2^r/„(a) cos2n^ 

n=l 

oo oo 

,/3sm2e _ j^(^) ^ 2 ^(-l)"/2„(/3) cos4r2^ ' ^Yl ^^""^^ sin2(2n + 1)^ 



= 27r/o(a)/o(/5) + 47r J](-l)"/2„(a)/2„(/5). 
An analytic continuation of the sum formula. 



n=l 



(31) 



n=l 71=1 

and integrating term by term. The cross-terms of the product of the two equations in Eq. |^ involve 
single powers of cos and sin and vanish upon integration. We are left with 



(32) 



CO 

JoiV + P'^ — 2a/3 cosy?) = Jo('^)-'o(/3) + 2 J„(q;) J„(/3) cosmp, 

n=l 



(33) 



at (y9 = 7r/2 simplifies Eq. |3^ to, 

ZV2 = 27r/o(0, ^ = v/^^+^. (34) 
Finally, the second integral in Eq. ^ can be computed by taking derivatives of Z^l'^: 

j (acos2^ + /3sin2^)e"'=°^2e+/3sin2e^^^ + ^i/2_ (35) 

Using these results, we arrive at the rotationally averaged energy Eejj given by Eq. ^ For large 
separation distances R, the effective interaction Ueff{R) = Eejf{R) — Eejf {00) defined in Eq. [Tl| can be 
expanded as in Eq. ^ where the coefficients are given by 

and 



A, 



2- A^^ 



d fhiO 



\IoiO 



d_ 



cos^ 2Q. (37) 
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Figure 1: Schematic of a protein inclusion. The top figure is a cut-away view of a membrane protein that 
contacts the continuum bilayer midplane on curve C. The contact angle on C is denoted 7 + ^7, while 
the bilayer deviation from a reference flat state is h{r). The bottom picture shows a possible ellipticity 
e in the projection of C onto the midplane. 





Figure 2: Two inclusions embedded in a local saddle deformation. The +/— correspond to 
raised/depressed regions of the membrane. The principle axis is aligned with the path joining the 
two raised regions (east-west). The principle axes of the inclusions {61,62) as well as the centerline 
joining their centers (fi) are measured with respect to this principle axis. 




Figure 3: (a). Rotationally averaged effective potential (Eq. [TsD as a function of protein separation in 
a flat membrane {Ht = Kh = 0) . (b). The minimum effective energy and its associated radius R*. The 
minimum of the potential is plotted as l/10|?7eif Note that R* quickly decreases when A increases 

above A* = ^2. For large A > 1, i?* ~ and \Ueff{R*) \ ~ 




Figure 4: Effective potentials between two inclusions embedded in anHj, — and constant Ki, membrane, 
(a). Ki, = 0.3; A = 0.1 for various fi. (b). Ki, = 0.2; A = 2.5, and (c). Kb = 2;A = 2.5. This latter case, 
although rare under physiological conditions, yields two energy minima which are physical manifestations 
of the qualitatively different minima depicted in (a), and (b). 




Figure 5: (a). The radii corresponding to interaction potential minima as a function of Ki, for A = 
0.5, 2, 4 and = 0, 7r/2. Curves that diverge signal a loss of the minimum (minimum radius R* — > oo) 
for parameters beyond those indicated, (b). The corresponding potential energy well depths at R* . The 
energies asssociated with A = 2; Q = and A = 2;fi = 7r/2 separate at K^, 1.1 when the Q = 
energy well disappears. The minimum energies associated with large A and Q = 7r/2 is still increasing 
for K, > 1.1. 



A-0.5;/<:,=0.1 
A=0.5;A:^=0.5 
A=2.5;ii:^=0.1 
A=2.5:ii:,=0.5 



A=0.5;/i:,=0.1 

A=0.5;/:^=0.5 

A=2.5;^:,=0.5 







71/8 



7C/4 



37U/8 



7U/2 



Figure 6: Angular dependence of (a). R, and (b). Ueff{R*,^) as functions of pair orientation angle 
Q. Minima arising mainly from background saddle (sensitive to Q) and ellipticity (insensitive to Jl) are 
shown. 
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Figure 7: (a). Second virial coefficient B2{A,Kii = 0,0.5,0.75,1). A negative virial coefficient is 
indicative of an overall attractive interaction such that the osmotic pressure is reduced from that expected 
in ideal solutions. The value 52(0,0) = tt^^'^/\/2 > corresponds to the virial coefficient of circular, 
repulsive {U — inclusions, (b). B2 as a function of background saddle for various ellipticity 

parameters A. 



